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Abstract 

We show soliton solutions of nonlinear Schrodinger equation on simple networks con- 
sisting of vertices and bonds, where the strength of cubic nonlinearity is different from 
bond to bond. We concentrate on reflectionless propagation of Zakharov-Shabat's solitons 
through a branched chain, namely, a primary star graph consisting of three semi-infinite 
bonds connected at a vertex. The conservation of the norm and the global current elu- 
cidates: (1) the solution on each bond is a part of the universal soliton solution on a 
simple 1-dimensional (1-d) chain but multiplied by the inverse of square root of bond- 
dependent nonlinearity; (2) nonlinearities at individual bonds around each vertex must 
satisfy a sum rule. Under these conditions, all other conservation rules for a simple 1-d 
chain have proved to hold for multi-soliton solutions on graphs. The argument is extended 
to other graphs, i.e., general star graphs, tree graphs, loop graphs and their combinations. 
Numerical evidence is also given on the reflectionless propagation of a soliton through a 
branched chain. 

1 Introduction 

Transport in networks with vertices and bonds [U [2] received a growing attention recently. The 
practical importance of this problem is caused by the fact that those networks mimic networks 
of nonlinear waveguides and and optical fibers [3] , Bose-Einstein condensates in optical lattices 
[I], superconducting ladders of Josephson junctions [HE], double helix of DNA j7], etc. 

Most studies so far, however, are restricted to solving the linear Schrodinger equation to 
obtain the energy spectra in closed networks and transmission probabilities for open networks 
with semi-infinite leads. 

On the other hand, with introduction of the nonlinearity to the time-dependent Schrodinger 
equation, the network provides a nice playground where one can see interesting soliton prop- 
agations and nonlinear dynamics in general. There already exist an accumulation of numer- 
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Figure 1: Primary star graph. Left panel: 3 semi-infinite chains connected at a vertex O; Right 
panel: space-discrete version of the left panel being used for numerical simulations in Appendix. 

ical studies of the soliton propagation through the discrete chain attached with small graphs 
[3 El QUI EH EE2] , where the discrete nonlinear Schrodinger equation(DNLSE) plays a role. How- 
ever, we see little exact analytical treatment of soliton propagation through networks, namely 
an assembly of continuum line segments connected at vertices, within a framework of nonlinear 
Schrodinger equation(NLSE)[13j. The subject is difficult due to the presence of vertices where 
the underlying chain should bifurcate or multi-furcate in general. A set of the continuity and 
smoothness conditions, which was exploited for the linear Schrodinger equation on graphs (see 
for example, Smilansky et al.[2]), would provide a natural boundary condition at each vertex. 
Then a soliton coming into the vertex along one of the bonds shows a complicated motion 
around the vertex. When time elapses, however, one can expect stable solitons with smaller 
amplitude propagating along each of bonds with one of them representing the reflected soliton 
along the incoming bond. The stable soliton in each bond is expected to be Zakharov-Shabat's 
solitons (ZSSs) discovered in the nineteen-seventies|14j . Hereafter the terminology of ZSS will 
cover from a single soliton through multi-soliton solutions presented in [14J. 

Interestingly, under an appropriate relationship among values of nonlinearity at individual 
bonds together with a suitable boundary condition at the vertex, one can see nonlinear dynamics 
of solitons with no reflection at the vertex. For example, let us consider an elementary branched 
chain (see Fig{T]), namely, a primary star graph (PSG) consisting of three semi- infinite bonds 
connected at the vertex O. 

As shown in FigfSJ our numerical simulation of DNLSE on the discrete version of PSG 
indicates: the soliton starting at lattice point x = 50 in the branch 1 enters the vertex at 
x = 200 and is smoothly split into a pair of smaller solitons in the branches 2 and 3, with 
neither reflection nor emergence of radiation at the vertex! Inspired by this discovery, we 
shall explore conditions among solitonic parameters and the relationship among strengths of 
nonlinearity at individual bonds in order to see the reflectionless propagation of solitons through 
networks or graphs. 

In this paper, we shall present the exact analytical treatment of soliton dynamics in net- 
works, by concentrating on the reflectionless propagation of ZSSs through vertices in the net- 
works. We shall search for the conditions for parameters characterizing the ZSSs and the sum 
rule for strengths of nonlinearity in each bond to satisfy the conservation rules for the norm 
and current during their propagation through networks. Network models we shall choose are 
star graphs, tree graphs, loop graphs and their combinations. We assume bonds of star graphs 
and edge bonds of tree graphs as semi-infinite, while other intermediate tree and loop bonds 
are taken finite. 

In Section II , using PSG, we shall show a basic idea of the soliton propagation along the 
branched chain. ZSS along the incoming chain is shown to bifurcate at the vertex, resulting in 
a pair of ZSSs with each propagating along the outgoing bonds. In Subsection II. A we establish 
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Figure 2: Time evolution of a soliton propagation through a vertex (numerical result): an 
example. Space distribution of wave function probability is depicted in every time interval 
T = 50.0 with time used commonly in branches 2 and 3. Abscissa represents discrete lattice 
coordinates defined in Fig{TJ Strength of nonlinearity at each bond are /3i — l,/3 2 = 1.5, /9 3 = 3. 
Initial profile is Zakharov-Shabat soliton in Eq.([H]) at t — with parameters a = 0.1, v = 0.1. 
Time difference in numerical iteration is At = 0.1. For the numerical method to solve nonlinear 
Schrodinger equation on primary star graph, see Appendix. 
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a connection formula for the current density at the vertex and find a relation among parame- 
ters characterizing the ZSSs. In Subsection II. B, we shall address the additional condition to 
guarantee the total current conservation rule, finding the sum rule for strengths of nonlinearity 
at each bond. In Subsection II. C, the boundary condition at the vertex is elucidated. In Sec- 
tion III, we show that the energy and all other conservations for a simple 1-dimensional (1-d) 
chain to hold for general solitons through PSG. In Section IV, we investigate the bifurcation of 
two-soliton and iV-soliton solutions at the vertex of a branched chain. In Section V, the cases 
of tree graphs, loop graphs and their combinations are investigated, where the soliton solution, 
sum rule and an infinite-number of conservation rules will be constructed by generalizing the 
result for PSG. Summary and discussion are devoted to Section VI. Appendix is devoted to the 
way of numerically solving the corresponding DNLSE on PSG. 



2 Single Soliton Propagation on Primary Star Graph 

2.1 Norm Conservation Rule and Connection Formula for Current 
Density 

We consider an elementary branched chain or PSG in the upper panel of FigH], where the 
vertex site is now taken as origin O. Space coordinates in individual bonds are here defined 
as bi ~ (—oo,0), b 2 ~ (0, +00) and b 2 ~ (0, +00). On each bond we have the nonlinear 
Schrodinger equation (NLSE) 

^ + ^ + A^| 2 ^ = 0, A: = 1,2, 3, (1) 

with Xk defined on — 00 < x\ < 0, < X2,x 3 < 00. It should be noted that the strength of 
nonlinearity (3k{> 0) may be different among bonds. The solution in PSG satisfies the following 
conditions at infinity: — > at x\ — > —00, ty k (xk) — > at x^ — > 00 for k = 2,3. One 
of the physically important conditions for the solution in PSG is the norm conservation. The 
norm is defined as 



N 



u 00 00 

|$|| 2 = J \ty 1 (x,t)\ 2 dx + J \^ 2 (x,t)\ 2 dx + J \^ 3 (x,t)\ 2 dx. (2) 



Let us find conditions for which the norm is conservative. For this purpose we calculate its 
t ime-der ivat i ve : 

d AT } d\^ x {x,t)\\ 7d\^,t)\\ , 7 d \M*M' 



—N = / 1 ^ ' " dx + 
dt J dt 



J t dx + S m dx ' < 3) 





From Eq. (CQ) we have the continuity equation, 



dt dx 
Using Eq.(jlj) in Eq.(j3j), we have 



dx 



(4) 



f f N = - 3l (0,t) + j 2 (0,t) + j3 (0,t), (5) 
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where j k (k = 1, 2, 3) is the current density defined by 



j k {x,t) = 21m 



%(x,t) 



dx 



From Eq.Q it follows that the norm is conservative only in the case, 

ji{0,t)=j 2 {0,t)+j 3 {0,t), 



(6) 



(7) 



which is the connection formula for the current density or the local current conservation con- 
dition at the vertex O. Similar conditions can be obtained for more complicated topologies. 
Below, Eq.flT]) will be evaluated explicitly. 

Let us assume that a single (bright) soliton in PSG is described with use of parts of ZSS 
([T3]) lying on individual bonds [k = 1,2,3) as follows: 



* k (x,t) 



a k 



V2 ex P 



if x 



sj~$ k cosh [a k (x + h~ v k t)] 



where v k , —l k and a k are arbitrary parameters characterizing velocity, initial center of mass and 
amplitude of a soliton, respectively. While these solutions are finite at the origin, it tends to 
zero at the infinity. 

We shall obtain the condition for the soliton solution Eq.([H]) to satisfy Eq.([7]). Noting 



dx 



«y - a k tanh(a fc (x + l k - v k t)) 



the current density becomes 

jk(x,t) = v k \^(x,t)\ 2 

Then Eq.Q is expressed as 
vial 



2v k al 



/3 fc cosh [a k (x + l k - v k t)\ 



v 2 a\ 



+ 



v 3 a\ 



/Sicosh [a^lx — vit)] /? 2 cosh [a 2 (/ 2 - v 2 t)] /3 3 cosh [a 3 (l 3 - v 3 t)\ 



(9) 



(10) 



Because of the linear independence of three functions, the connection formula Eq. ffTTT) is 
satisfied only in the following three cases: 



m 



[i) a 3 v 3 = 0, ai/i = a 2 / 2 , a-i^i = 

ii) a 2 v 2 = 0, aih = a 3 / 3 , a\V\ = 
ai/i = a 2 l 2 = a 3 / 3 , a x v x = a 2 v 2 



ai a 2 

ai a 3 
ai a 2 a 3 



(12) 

(13) 
(14) 



The above results are also available more straightforwardly by demanding the norm Eq.([2]) to 
be time- independent. The cases (i) and (ii) suggest that a soliton simply moves from the bond 
b\ to bonds b 2 and 6 3 , respectively, while the case (iii) is indicative of the splitting of a soliton 
on b\ into two parts with one appearing on b 2 and the other on 6 3 . Besides Eqs.( fT2|) -( fT4T) . there 
are additional constraints on solitonic parameters and strength of nonlinearity by noting the 
global current conservation, which will be described below. 
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2.2 Current Conservation Rule and Sum Rule for Strength of Non- 
linearity 

We then demand the total current of the soliton in PSG, which can be calculated as 

+00 +00 

J= J ji(x,t)dx+ J Mx,t)ds+ J js(x,t)d X (15) 
-00 

= E(-i) 2fc-1+lH ir tebWk - v k t)) + j2 2 ^p. (is) 

k=i i Jk k=i P k 

It follows from Eq. (fl6l) that the total current is conserved, if the additional conditions 

ViOi _ v 2 a2 n? x 

Pi ~ P2 ' 1 } 

v\ai v 3 a 3 

IT "IT' (18) 

and 

VlOi _ V2CI2 ^3«3 / 1Q s 

Pi " & A ' 1 j 

are satisfied for cases (i), (ii) and (iii), respectively. 

The norm conservation rule in Eqs. (fT2]) - (fT4j) and current conservation rule in Eqs. f)17l) - ffT9l) 
can be simultaneously satisfied under the following cases: 

[i] a 3 = 0, a x = a 2 = a; vx = v 2 = v; h = k = I, Pi = P2 = 0, (20) 

[ii] a 2 = 0, ax = a 3 = a; Vx = v 3 = v; h = h = I, Pi = fh = P, (21) 
[iii] ax = a 2 = a 3 = a; v t = v 2 = v 3 = v; h = k = h = h (22) 

irk + h (23) 

where a,v, I and P are arbitrary constants. Equations (!20l) - (l22|) imply: Firstly, the soliton solu- 
tion at each bond should be a part of the identical ZSS, namely the main solitonic parameters 
should be common to individual bonds (line segments) in the graph except for the strength 
of nonlinearity. In particular, solitons on b 2 and b 3 are initially located outside these bonds; 
Secondly, the soliton can bifurcate in passing through the vertex if the strengths of nonlinearity 
Pj at individual bonds bj satisfy the sum rule in Eq. (1221) ! This rule can also be obtained in a 
different way with use of symmetry argument in solving DNLSE (see Appendix). 

To be explicit, we shall see the following dynamics: In the first two cases, [i] and [ii], the 
soliton coming first from the bond bx disappears or becomes a ghost at time r = -, when a new 
soliton appears in either one of b 2 and b 3 . In these cases, we have Px = P 2 = P or Px = P 3 = P, 
and therefore the soliton propagation is nothing but that in an ideal 1-d chain. 

The third case [iii] is the most interesting, where the soliton at bond bx splits into two parts 
and appears in both of 62 and b 3 , as shown in Figj3l This is a novel feature of the soliton 
propagation through a branched chain and networks in general. Precisely speaking, the soliton 
dynamics here is governed by a single characteristic time r = -. While for < t < r the soliton 
at bx is a real one and those at 62 and 63 are ghosts, for r < t the soliton at 61 is a ghost and 
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t = t = t, 

Figure 3: Splitting of Soliton. ti > r. Broken curves represent ghost solitons. 

those at 62 and 63 are real. Common to the cases [i]— [iii] , the incoming real soliton on b\ and 
outgoing ghost solitons at 62 and b% arrive at the vertex O at time r simultaneously. 

In a similar way, the time-reversal process is obvious that two solitons start to move from 
bonds 62 and 63, meet each other at the vertex at time r and move as a single soliton along the 
bond b\ towards —00. 

2.3 Boundary Condition at Vertex 

In closing this Section, we should clarify the nature of the boundary condition at the vertex. 
Equations (!20l) - (l22l) have led to the issue that the solution on each bond is a part of the 
universal soliton solution on a simple 1-d chain but multiplied by the inverse of square root of 
bond-dependent nonlinearity. Indeed, under the condition in Eq. (1221) . the soliton solution in 
Eq. ([8]) can be written as 

Vk(x k ,t) = ^-iq(x k ,t), (24) 

where the functions q(xi,t) and q(x2 t 3,t) are restricted to (— 00; 0] and [0, +00), respectively 
and stand for the individual part of the /5-independent universal solution q(x,t) which satisfies 
the nonlinear Schrodinger equation with (3 = 2, 

iqt + Qxx + 2q\q\ 2 = 0, —00 < x < +00. (25) 

Here, *(#, t) itself is neither continuous nor smooth at the vertex. As for q(x), however , there 
is no singularity there: lim a . 1 _»_o q^ n \xi) = lim^ 3 ^ +0 q^(x 2t 3) for any n = 0, 1, 2, • ■ • with q^ 

1 /2 

the n-th order derivative of q. Equivalently, the solution scaled by by j3 k (A: = 1,2,3) satisfies 
the boundary condition at the vertex: 

lin n /f*S B) W= hm ^^g^s), (26) 

21— >— U 22,3^+U 

for any n = 0, 1, 2, • ■ • with ^ the n-th order derivative of ^k- The absence of singularity 
of any kind in the scaled function f3^ {%k) is the reason why we see neither reflection nor 
emergence of radiation at the vertex, as is evidenced numerically in FigfJ (see also Appendix). 

This implies that the boundary conditions used in the present work are different from those 
in the case of time-dependent linear Schodinger equation on graphs (for example see, Kottos & 
Smilansky 1997). If we would take ordinary continuity and smoothness conditions, e.g., 

**(0,*) = *i(0,t), (27) 
9*1 (9*9 9*3 

^ CM =^0,* ) + 0,f , (28) 
ox ox ox 

we shall see a completely different nonlinear dynamics of solitons such as reflection of a soliton 

at the vertex. The initial value problem under the conditions in Eq. (l28|) at the vertex will be 

treated elsewhere. 
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3 Energy and Other Conservation Rules 

We now proceed to the calculation of the energy of soliton in the graph. The total energy can 
be evaluated as 



E 



0*i (M) 



dx 
d^ 2 {x,t) 



dx 



+00 / 

/ 




\ 


dx 



^ 1 ( Xi t)\ 4 \dx + 



^\* 2 (x,t)\ 4 \dx 4 



-^|^ 3 (x,t)| 4 I dx. 



(29) 
(30) 
(31) 



With use of Eqs.flED an d ©> we have 



dx 



-%\M*,t)\'=(4-*i\ + 



2a\ tanh 2 [afc(x + l k — v k t)\ 



cosh [a k (x + l k - v k t)} 



(32) 



(33) 



In the most essential case [hi], for example, Eqs. (l3"Tj) and fl3"3"|) together with Eq.(j2"2"|) lead to: 



E = 2a 



v 2 \ 2a 2 
! — — a 2 J tanh(a/ — avt) H — — tanh 3 (a/ — avt) 



(34) 
(35) 



which proves constant under the sum rule in Eq.(l22~l). 

One can generalize the argument so far beyond a single-soliton solution: So long as the 
general solution on PSG is described by parts of the corresponding universal scaled function 
q(x, t) as shown in Eq. (1241) . all the conservation laws for 1-d chain hold for PSG under the sum 
rule Eq. (l22l) . Applying Zakharov-Shabat's theorem ([14J), we obtain the general conservation 
rules 



(2i) n C n p^ = E A-' / fn(q k (x k ,t))dx k , 



where C n are constant, and f n (q(x)) obeys the recursion relation (see Eq.(35) of [14"]): 



fn+l = (-fn) + E fjfh fl = W 

\" / j+l=n 



(36) 



(37) 



In fact, with use of Eq. (|22|) . the r.h.s. of Eq. (1361) turns out: 



Pi 1 I f n (q(x,t))dx+((3 2 - 1 + /3^) J f n (q(x,t))dx 



+00 

Pi 1 I fn(q(x,t))dx = fc\2i) n C n , 



(38) 



where the second equality comes from the conservation rule for the 1-d chain (|14j). 

It is easy to see that f n is the 2n-th order polynomial of q and its derivatives with respect 
to x, written in the following form 



fn ^ ^ bs-Pn,2s(qi Qxj Qxxi 



(39) 



s=l 



where P n ^ s = q kl (q*) k2 q k3 (q*) k4 • ■ • with k% + k 2 + k 3 H = 2s. 

In this way one can obtain an infinite number of conservation laws in PSG, 



(2i) n C n ^ 1 = fracU E / E M y ) ^n,2 S (* fc , • • -)dx k 



(40) 



k=1 b k s=1 



In Eq.f l40p . the cases n = 1,2, and 3 give the norm, current and energy conservation rules 
in Eqs.(T5]), (ITSj) and (13~T1) . respectively. Some higher-order conservation rules are as follows: 



k b k 



d 3 ^ 

dx\ 



dx k 



^k\ 2 (x k ,t)dx k , 



(41) 



(2i) 5 C B ^ 1 



-E 



-3ft 



a 2 ^ 



+ $i*.|6_& 



_d_ 

dx k 



5xa 



(x k ,t)dx k . 



(42) 



(43) 



The above treatment is also true for more general star graphs consisting of N semi-infinite 
bonds connected at a single vertex. In such cases, the initial soliton at a bond splits into N — 1 
solitons in the remaining bonds, and the extended version of Eq. (122]) is: 



I (7 = 1,2,. ..,JV), 7T=Ey 



= a, fj = u 3 lj 



(44) 



where a, u and I are arbitrary constants. 

All the soliton solutions by Zakharov and Shabat, as they stand, are applicable to networks, 
but, in order to see the bifurcation of the soliton solution at vertices, the strengths of nonlin- 
earity at individual bonds should be different and satisfy the sum rule Eqs.( |22i) or ( H4"l) at each 
of vertex, which we shall see in details in the following Sections. 



4 Mult i- Soliton Solutions 

For the case when the condition like Eq. fl2"2l will be satisfied, one can construct more general 
multi-soliton solutions of NLSE satisfying all conservation rules. Below we shall give such 
multi-soliton solutions on PSG. 
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4.1 Soliton Collision Described by Double-Soliton Solution 



Here, three solitons incoming along each of 3 bonds towards the vertex are shown to be scattered 
to three directions. The following solution of NLSE on PSG describes double-soliton solution 
or a soliton collision: 



= iq(xx,t), xiG(-oo,0], 
Pi 



-iq(x k ,t), x k e[0,+oo), k = l,2. 



(45) 
(46) 



Here 
where 

with 



q(x, t) = A+(x, t) exp(— iz\ x) + A-(x, t) exp(—iz*ix) 



A 



± 



D ± B TT - D T B ±T 
B++B — — B j — B — l 



D± = — c*± exp(— iz* ± x) 



v Cl cj 



— exp(i(2zi - z\ - z*L)), 



B 



±± 



E 



— exp(z(2^-24))-l, 



z=+ ,_ {zi - Z±j 

z± = ^± + if]± = const, 
c± = c±(0) exp(4iz±t), 
where £±,r)± and c±(0) are scattering data characterizing the solitonic parameters. 





t = o r = ^ 

Figure 4: Collision of solitons. £i > r. Broken curves represent ghost solitons. 

A collision among solitons is possible for < 0. It follows from the above double- 

soliton solution that after a collision two solitons penetrate each other and are split into two 
independent solitons. Velocities of solitons are characterized by the parameters £±, while its 
amplitude depends on the parameters, /3&, and n±. On the bond b%, for example, the double 
solitons are regarded as consisting from a real one with velocity v + > on the bond b\ and a 
ghost one with velocity V- < on the ghost bond extended from b\ to the positive x region. 
Let consider the case that a pair of real and ghost solitons meet at the vertex at time r, by 
properly choosing the initial position of the centers of mass which depends on the scattering 
data. Then we can see collision among three real solitons in FigJH While for t < r the real 
soliton on b\ (^2,3) has the velocity v = v + (v = for t > r it acquires v = w_ ( v = v + ), 
describing a collision among 3 real solitons. 
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4.2 iV-Soliton Solution. 

Similarly to the case of two-soliton problem we can obtain iV-soliton solution of the NLSE on 
PSG by assuming Eqs.( l22l) and (1241) . Let q(x,t) be the iV-soliton solution of NLSE with the 
nonlinearity (3 = 2 in Eq. (1251) . Then the iV-soliton solution of Eq.([T]) on PSG can be constructed 
by multiplying the corresponding universal solution g(x, t) with the bond- dependent factor y-J^i 
on individual bonds b^. 




Figure 5: Collision of multi-solitons, which apparently breaks conservation of particle numbers. 
t 1 >r. N = 3,M = 2. 

It should be noted that at arbitrary moment of time the bonds 62 and 63 have the same 
number of solitons, M. However, at the same time the bond b\ has N — M solitons, as 
exemplified in FigJSJ With use of iV-soliton solutions, a variety of splitting of solitons at the 
vertex is found depending on the initial velocities, which apparently breaks the conservation of 
particle numbers. As already described in Section IV, all the conservation laws hold for this 
solution. 

5 Other Types of Graphs 

Now we proceed to explore soliton solutions of NLSE on other kind of graphs and explore the 
sum rule and conservation rules for solitons to propagate through these graphs. 




Figure 6: Tree graph. b\ ~ (— oo,0),6n ~ (0,L n ),f> 12 ~ (0,L 12 ), and buj ~ (0, +00) with 

l,J = l,2,--: 

An example of the graph for which the soliton solution of NLSE can be obtained analytically 
is a tree graph in Fig. [61 Hereafter, for an arbitrary one of bonds in the tree graph, we shall 
employ an abbreviation like br = buj.... On each bond br we have NLSE given by Eq. ([!]) and 
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for each vertex the following conditions is satisfied: 

T = ^h (47) 

which is again available from the norm and current conservation rules. The soliton solution 
satisfying these conditions can be written as 

ty r (x r ,t) = ^J^-iq(x r + s r ,t;s r ), x r eb r . (48) 

Here paremeter sr is the length of the path that soliton passes from b\ through br- For tree 
graphs this parameter is given as 

Si = sii = I, Suj = I + Lu, 



s r = suj-im — I + Lu + L U j + ■ ■ ■ + Luj...i, (49) 

where —sr represents an initial location of the solution (arbitrary part), and Lu, Luj, ■ ■ ■ , Luj...i 
are lengths of the finite bonds prior to br = buj-im- 

Below, applying the induction method, we give a proof of conservation rules for soliton 
solutions of NLSE on any tree graph. Let us denote the tree graph in Figj6] as G and assume 
the conservation rules to hold in G: J2b r eGflr I fn{q{ x v + sr,t))dxr = (2i) n C n p 1 ' . Then we 

construct an enlarged tree graph in the following way: First, cut an arbitrary one of the right- 
most semi-infinite chain 6a ~ (0, +oo) at a point A located by distance La from the nearest 
vertex and then attach iV semi-infinite bonds to the point A which now becomes a new vertex 
point. Namely the bond 6a is now replaced by the finite bond 6a ~ (0, La) connected with iV 
semi-infinite bonds b\ m ~ (0, +oo) with m = 1, • • • , N . The enlarged tree graph thus obtained 
is denoted as G'. In the same way as in Eq.( l36|) . the general conserved quantity for G' is given 
by 

E breG _ 6A ^r 1 / fn(q(x r + s r , t))dx T + p£ x J fn(q(x A + s A , t))da 

br 'b A 
N . 

,-i — 1 In// t ,\\ 

IX An 



IX \ 



+ E ^Am / fn{q{Xhm + S\ t))d 
m=l ~ 

= E breG ^r 1 / fn(q(x r + s r , t))dx T 

b r 

+oo N +oo 

-Pa 1 J fn(q(x + s A , t))dx + J2 Pa™ J fn(q(x + s A , t))d 



m=l 



V m =l / £ 



f n {q{x + s A ,t))dx. 



Here J2b r eG-b A an d J2b r &G imply summations over all bonds in G except for 6a and over all 
bonds of G, respectively. It is clear that the final expression becomes constant (2i)"C n/ 9 ] " 1 
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Figure 7: Graph with loops. b ~ (— oo, 0), b n+ i ~ (0, +oo), ~ (0, L) with k — 1, 2, n. 

under the sum rule in Eq. (H7|) . Thus, starting from PSG in Fig. [TJ and repeating the above 
procedure, we can get the conservation rules for all possible tree graphs. 

Another example for which soliton can be easily obtained is a graph with loops (see FigJTJ). 
This graph consists of two semi-infinite bonds whose edges are connected with n bonds having 
finite lengths. Again, requiring the following conditions for the coefficients of NLSE: 

1 = " 1 = 1 

we can write the soliton solution by Eqs. (l4*8l) . 

Also, the exact soliton solution can be obtained for the graph in FigJH] where the corre- 
sponding condition for the parameters, (3^ is required. This graph can be considered as a loop 
graph connected with 3 semi-infinite bonds. 




b 6 



Figure 8: Loop with semi-infinite bonds. &i ~ (— oo, 0), 62, 63 ~ (0, +00), bk ~ (0,Lfc) with 
fc = 4,5,6. L 6 = L4 + L5. 




g 2 



Figure 9: Combination (right panel) of loop G\ and tree G2 graphs (left panel). 

Combining the above topologies one can construct different graphs(having semi-infinite 
bonds) for which soliton solution of NLSE can be constructed. To do this consider two graphs G\ 
and G2 where there exist soliton solutions given by Eq. (1481) . Soliton solution on the combination 
of the graphs G\ and G 2 can be constructed using one of two methods below. 
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Method I. Let A be a point on a finite bond of the graph G\ and B be a point on a semi- 
infinite bond of G 2 . Connecting two graphs by putting together the points A and B, we reach 
the graph in Fig. 



Method II. Fixing two points Ai and A 2 on a bond of graph G\ and cutting the part between 
these points and doing the same thing with the points B\ and B 2 of bond in the graphs G 2 we 
can connect two graphs by putting together the points A\ , B\ and A 2 and B 2 , respectively, 
and thereby can reach a new graph in Fig. [TUJ 

In these different types of graphs, the soliton solution of NLSE can be constructed under 
the conditions given by Eqs. (l4"Tl) and ( 1481) . It should be noted that throughout in our approach 
the graphs are supposed to have at least two semi-infinite bonds. 

6 Summary and discussions 

We have explored soliton solutions of nonlinear Schrodinger equation (NLSE) on simple net- 
works. We first concentrated on reflectionless propagation of a Zakharov-Shabat's soliton 
through a branched chain, namely, a primary star graph. To satisfy the conservation of the 
norm and global current in the graph, the solution on each bond should be a part of the 
universal soliton solution on a simple 1-d chain but multiplied by the inverse of square root 
of bond- dependent nonlinearity. Besides this, nonlinearities at individual bonds around each 
vertex must satisfy a sum rule: the inverse nonlinearity at an incoming bond should be equal 
to the sum of inverse nonlinearities at the remaining bonds. Under these conditions, all other 
conservation rules for solitons in a simple 1-d chain have proved to hold for solitons propagating 
through graphs. With use of Zakharov-Shabat's two-soliton solutions we also find a collision 
among three solitons at the vertex, and with use of iV-soliton solutions, a variety of splitting 
of solitons at the vertex was found depending on the initial velocities, which apparently breaks 
the conservation of particle numbers. The argument is extended to general star graphs, tree 
graphs, loop graphs and their combinations. To see all conservation rules to hold, a set of in- 
verse nonlinearities should satisfy the generalized sum rule at each of vertices, which we proved 
by the induction method. 

So long as the sum rule for strength of nonlinearity at each vertex holds, the boundary 
conditions (: connection formulas) there for scaled wave functions are normal, and there is no 
singularity at vertices that generates reflection and radiation. Although there exist accumula- 
tion of analytical studies on initial value problems on the semi-infinite chain (Ablowitz & Segur 
1975; Fokas et al. 2005) and the finite chain (Ramos & Villatoro 1994; Fokas & Its 2004), no 




Figure 10: Same as FiglU but a different combination. 
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corresponding ones in networks or graphs have appeared up to now. Under the boundary con- 
nections different from the present paper, soliton dynamics would become more complicated, 
whose analysis is a next challenging subject. 
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Appendix. Numerical method to solve nonlinear 
Schrodinger equation on primary star graph 



With use of space discretization 



,d 2 jj 
dx 2 



ipi+\ — 2'0 i + '0j_ 1 ), the nonlinear Schrodinger equation 



in the 1-d continuum with neither branches nor vertex 

.difj d 2 ijj 



dt 



dx 2 



can be reduced to 



Mi 
' dt 



+ (2 - i3\ipi\ 2 )ipi - Vi+i, 



(50) 



which is rewritten in a matrix form as 




''dt 



tpi-i 
V : / 



H(t) 



( \ \ 
V i / 



(51) 



with 



H(t) = 





w 




-1 2-/3|^_i| 2 


-1 





0-12 


-m\ 2 


-1 





-1 2 





V 



(52) 



Then, by carrying out the time discretization with time difference At, Eq. (j51j) reduces to 



^x(t + At) 

Mt + At ) 
ijj i+1 {t + At) 



V 



exp(-iH(t) ■ At) 



( \ \ 



) 



\ 



(53) 



/ 
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It is obvious that Eq. (!53|) conserves the norm because of the unitarity of exp(—iH(t) ■ At). It 
is sraight-forward that the diagonalization 



p- x H{t)P 



(e x • • ■ 

e 2 



gives rise to 



exp(-iH(t) ■ At) = 

( exp(-ze x At) 

exp(— ie 2 At) 







V 



P- 1 



(54) 

(55) 
(56) 



In case of a branched chain, i.e., a primary star graph (PSG), we consider its discretized 
counterpart and introduce the numbering as in FigJTJ Equations to generalize Eq. (1501) are given 
by 



(1) 



dt 



-^ ) 1 + (2-/? 1 |# ) l 2 )# ) -^ 



+1' 



E) / ( 2 ) 

^--^ 1 + (2-/? 2 |^ ) | 2 )^ 2) -^ (2 + ) 



dt 
dt 



-^ + (2-/?3|^ 3) | 2 )4 3) 



(57) 
(58) 
(59) 



which correspond to bonds bi, b 2 and b 3 , respectively. 

The important problem is to search for the connection formula at the vertex, which will be 
resolved as follows. Let call the end of b± as K site. Similarly the starts of b 2 and 63 are taken as 
L and M sites, respectively. Introducing virtual wave functions Vx-i an d vm-i an d establish 
for their relationship with . As a manifold of the global solution on PSG, we assume a 



discretized version of Eq.( l24j) : 



VWk 



ig{xj,t), 



(60) 



where k(= 1, 2, 3) denotes individual bonds and the discrete lattice variable i runs over PSG in 
Fig{TJ Because of the continuity of g(xj,t) at the vertex, we obtain a connection formula: 



(61) 



On the other hand, with use of suitable parameters s 2 and S3, a virtual wave function ip^+i 
should be 



1 



S2 + S3 



S2\ 







Then Eqs.fl57l)-(l59~l) at the vertex can be explicitly rewritten as 



.dip 



(i) 

K 



dt 



-^ ) _ 1 + (2-/3 1 |^ ) |>«- 



s 2 i 



' dt 



S2 + S3 

£*P + (2-Ah>£?|')*2>-*2> 
Pi 



M+l' 



(62) 

(63) 
(64) 
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Lining up ^\ i/jf^ and ifj^ vertically and rewriting Eqs. (l571) - (l59l) with Eqs. (l63l) - (1651) in 
a matrix form, we obtain the equation like Eq. (15 11) with Eq. (1521) . but with a modified real 
matrix H(t). To conserve the norm, H(t) should be symmetric, which imposes the following 
relationship: 

S2 [fa _ fjh s 3 //% _ /ft 

S2 + S 3 V A V ^2 ' S 2 + S 3 V Pi V /V 

As a result, we have the sum rule for three kind of strengh of nonlinear ity: 

irk + h (67) 

which agrees with Eq.( l22l) obtained from the norm and current conservation rules for the PSG 
in the text. Using Eq. (l66|) , Eq. (l63!) can be replaced by 



(i) 



^ + (2 - - (M*%> + M^s) ■ (« 



dt lift Ift 



3 / 

By numerically solving Eqs. fl571) - fl59l) with Eqs. fl68|) . fl6%|) and fl65|) under any initial condition, 
one obtains nonlinear dynamics of solitons without reflection at the vertex. Figure [2] is obtained 
under the initial profile in Eq.Q with Eqs.(|67l) or ( |22l) . 
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